Quantum Transport from the Perspective of Quantum Open Systems 
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By viewing the non-equilibrium transport setup as a quantum open system, we propose a reduced- 
density-matrix based quantum transport formalism. At the level of self-consistent Born approxi- 
mation, it can precisely account for the correlation between tunneling and the system internal 
many-body interaction, leading to certain novel behavior such as the non-equilibrium Kondo effect. 
It also opens a new way to construct time-dependent density functional theory for transport through 
large-scale complex systems. 
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Conventionally, quantum transport and quantum dis- 
sipation are treated with different approaches. For in- 
stance, the former (in mesoscopic context) is usually de- 
scribed by the Landauer-Biittiker theory and the non- 
equilibrium Green's function (nGF) approach 0, 0, 
whereas the latter by the reduced density matrix equa- 
tion 0- Nevertheless, both are quantum open systems, 
with cither the non-equilibrium electron reservoirs (elec- 
trodes) or the dissipative thermal bath as their environ- 
ments, as schematically shown in Fig. 1. ft is thus of great 
interest to develop a unified language to bridge them. 

This motivation can be historically dated back to 
the phenomenological rate equation and quantum Bloch 
equation approaches to transport 0, There, either 
implicitly or explicitly, the electrodes are viewed as dissi- 
pative reservoirs. Along this line and based on our work 
in quantum measurement of solid-state qubit 0, 0] , we 
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FIG. 1: A unified picture for (A) quantum transport and (B) 
quantum dissipative open system, where the transport system 
is regarded as the "system of interest" , and the electrodes as 
"environmental baths" . 



developed recently a quantum master equation approach 
to quantum transport |§j]. The reduced dynamics in- 
volved there was originally constructed under the cumu- 
lant second-order expansion (Born approximation). In 
this letter, we re-formalize it in the spirit of self-consistent 
Born approximation (SCBA), in order to make the for- 
malism not only convenient but also accurate enough 
in practice. Moreover, by reducing the many-particle 
density matrix formalism to single-particle one, we will 
also construct an efficient approach for large-scale (e.g. 
molecular) system applications. 

The typical transport setup, see Fig. 1(A), can be de- 
scribed by 
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Hg is the system Hamiltonian, which can be rather gen- 
eral (e.g. including many-body interaction), aj^ (a M ) is 
the creation (annihilation) operator of electron in state 
labelled by "/x" , which labels both the multi-orbital and 
distinct spin states of the transport system. The sec- 
ond and third terms describe, respectively, the two (left 
and right) electrodes and the tunneling between the elec- 
trodes and the system. 

To contact with the quantum dissipation theory for 
quantum open systems, let us introduce the reservoir 
operators = ^ak^kda^k = Jl m + fn^- Ac- 
cordingly, the tunneling Hamiltonian H' reads H' = 
(ajj-F/i + H.c.) . Treating H' perturbatively up to the 
cumulant second-order expansion, it yields [tj 
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p{t) = -ic p (t) - / dT{£{t)g{t,T)c'(T)g^{t,T))p{t). 

Jo 

The reduced density matrix is defined as p(t) — 
Tre[pT(t)], by tracing out the reservoir degrees of free- 
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FIG. 2: Diagrammatic illustration for the second-order tun- 
neling self-energy processes, on the Keldysh contour. The up- 
per and lower horizontal lines describe the forward and back- 
ward propagation of the transport system, which is treated 
exactly in terms of the system Green's function G(t,r). The 
dashed line stands for the tunneling between the system and 
electrodes. 



t, can be carried out via dt'a„(t,t') = iS(t — t')a v + 
i[Hs(t'), a v (t, t')]. Thus, the time integral in Ap , which 
becomes now the type of ~ J * dt'C&pl(t — t')a v (t, t'), can 

be calculated accordingly. Inserting the obtained 
into Eq. J3J, the time-dependent phenomena associated 
with either quantum dissipative dynamics or transport 
current can be easily treated. For clarity, we hereafter 
assume the system Hamiltonian to be time- independent, 
unless further specification. 

Now we consider the possibility to go beyond the 
second-order self-energy process diagrammatically shown 
in Fig. 1. An efficient scheme follows the idea of the 
well-known self-consistent Born approximation (SCBA), 
i.e., the free propagator defined above, IB™ (i) = 
—iQ(t)Q(t, 0), is replaced by an effective propagator H(t). 
The latter is obtained formally via the Dyson equation 

M 



ri(i) = -is(t) - icn(t) 



- i')n(t'), (4) 



dom from the total system-plus-reservoirs density matrix. 
The Liouvillian superoperators are defined as C{- ■ • ) = 
[H S , (•••)], £'(•••) = [H', (■■■)}, and £?(t,r)(-) = 
Git, r)(- • • )G>(t,T) with G(t,r) the usual propagator 
(Green's function) associated with the system Hamilto- 
nian H$- 

The integral kernel in Eq. @ , which is in the so-called 
artial ordering prescription (POP) (or time-local) form 
describes the second-order tunneling self-energy. At 
the second-order level, one may replace p(t) in the last 
term of Eq. Q with Q(t, r)yo(r), leading the tunneling 
integral kernel to (£'(t)£(t, t)£'(t))p(t), being in the 
chronological ordering prescription (COP) (or memory) 
form 0. The corresponding four terms in the conven- 
tional Hilbert space 01, depicted on the real-time 
Keldysh contour in Fig. 2, provide a clear diagrammatic 
interpretation for the second-order tunneling self-energy 
process. 

Explicitly tracing out the states of electrodes, Eq. @ 
gives rise to 
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For time-independent system Hamiltonian, A^ = 
Ea=i l «^ ) = J2 av JZo dtC ^(Tt)[iau{t)l with 
C { a +l(t) = </t M (i)/ Q „(0)), CL-I(t) = (f a ,(t)fU0)), and 
a v (t) = -i@(t)g(t,0)a u = n(°'(t,0K. Note that the 
step function O(i) extends the lower bound of the time 
integral from to — oo, whereas the extension of the up- 
per bound from t to oo results from the consideration of 
Markovian approximation. 

For time- dependent system Hamiltonian, the timc- 
translational invariance breaks down, we thus define 
a„(i,i') = — i®(t — t')Q(t,t')a„. The backward evolu- 
tion of a,u(t,t') with respect to t', starting from t' = 



or n(w) = [lo — C — 'E(lo)]^ 1 in frequency domain, with 
S being the irreducible self-energy defined again by Fig. 
2. Accordingly, fi„(i) = U(t)a v , and 
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Equations 0) constitute the basic ingredients of the 
proposed SCBA scheme. This type of self-consistently 
partial summation correction has included an infinite 
number of higher order tunneling processes into the re- 
duced system dynamics. The resulting non-trivial effect 
on quantum transport will be demonstrated soon. 

So far, the trace is performed over all the electrode 
states, and the resulting Eq. J3J) is a unconditional master 
equation for the system reduced dynamics. To character- 
ize the transport problem, we should keep track of the 
record of electron numbers entering the right reservoir 
(electrode). Following Refs. and|a one can obtain a 
conditional master equation for the reduced system den- 
sity matrix, p^(t), under the condition that n electrons 
have arrived at the right electrode until time t. On the 
basis of p( n \t), one is readily able to compute various 
transport properties, such as the transport current, the 
probability distribution function P(n,t) = Tr[p( n )(t)], 
and the noise spectrum 0. For transport current, it can 
be carried out via I(t) = eJ2 n nTr[p( n ) (t), giving rise to 



I(t) 
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p(t)+R.c. . (6) 



Compared to other transport formalisms, Eqs. ©-© 
provide a convenient framework for quantum transport. 

As an illustrative application, we consider the non- 
trivial problem of quantum transport through strongly 
interacting quantum dot, under the well-known Ander- 
son impurity model Hamiltonian: Hs = S/ t ( e o a I a M + 
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^n^np). Here the index fx labels the spin up ("t") and 
spin down ("J") states, and p, stands for the opposite spin 



orientation. The electron number operator n M 
and the Hubbard term Un^n^ describe the charging ef- 
fect. Apparently, the reservoir correlation function is di- 
agonal with respect to the spin indices, i.e., C^i(t) = 
5pvC£l(t), and C&l(t) = £ fe \ t^e+^n^ (e k ). 
Here (e&) = n a (efc) is the Fermi distribution func- 
tion, and ria \ek) = 1 — n a [ek)- Accordingly, we 

have Ahtp = T afl J |^rtL ±) (e)[in(e)]a / _ I , where, under the 
wide-band approximation, we have introduced r QM = 
2Trg a \t afJ- k\ 2 , and assumed it energy independent. From 
Eqs. © and JjJJl, the stationary current is obtained as 

I = J ^Im[H(6)][n L ( e ) - n R {e)]. (7) 

For the single level system under study, the propagator 
in energy space simply reads n(e) = [e — eo — S(e)] -1 . 

Within the SCBA scheme, the self-energy E can be 
explicitly carried out via Fig. 2. However, in the case 
of strong Coulomb repulsion, the dot can be occupied 
at most by one electron. As a result, it can be eas- 
ily proven that only Fig. 2(C) and (D) contribute to 
the self-energy. Physically, replacing the bare system 
propagator with the effective propagator corresponds to 
including the infinitely multiple forward and backward 
tunnclings between the system and the same electrode. 
This is in fact a tunneling-induced quantum fluctuation, 
which would lead to the level broadening and the non- 
trivial interference between tunneling and system inter- 
nal interaction. Explicitly, in large-?/ limit, the real 
and imaginary parts of the self-energy read ReE(e) 
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tively [ill Il2|. Here = \/{kgT) is the inverse tem- 
perature, fj, a the chemical potential of the electrode, tp 
the digamma function, and to denotes the spin degen- 
eracy, (i) For to = 1, i.e. neglecting the spin degree of 
freedom, ImII(e) gives the well-known Breit-Wigner for- 
mula, which appropriately includes the level broadening 
effect, (ii) For to > 2 (e.g., to = 2 for spin 1/2), the above 
self-energy correction would result in rich behaviors, de- 
pending on the relative values of the parameters such as 
the temperature and the position of eo with respect to 
the Fermi levels. Detailed discussions, in particular the 
non-equilibrium Kondo effect, are referred to literature, 
e.g. Refs.mil 

Application to Large Scale Systems. — By far, the 
transport-related density matrix formalism has been con- 
structed in many-particle Hilbert space, which may re- 
strict its direct application only in small systems. For 
large-scale systems in the absence of many-electron in- 
teraction, we first recast the formalism to a very sim- 
ple version in terms of the reduced single-particle density 
matrix (RSPDM), a^t) = Tr[ata M p(t)], which greatly 



reduces the dimension of Hilbert space, thus saves com- 
puting expense. To account for the electron-electron in- 
teraction, we then propose an efficient time-dependent 
density functional theory (TDDFT) scheme. Note that 
it is quite natural to combine the TDDFT technique with 
the present RSPDM formalism, since the former self- 
consistently amounts to the many-body interaction but 
still keeps the single-particle picture |l4j . 

(i) Time Independent System Hamiltonian: For 
simplicity, we first proceed our derivation in the 
single-particle eigenstate basis, which is denoted as 
• • • }• In this representation, Aajx — 
£„ C^l^e^a^, and the equation of motion for the 
RSPDM can be readily obtained by applying Eq . di- 
rectly for o-^it) = Tr[a£a M/ o(t)]. We have fTMlll| 



-i[h,a]-\{[c { 



H 



(8) 



Here, h is the single-particle Hamiltonian or the Fock ma- 
trix within the TDDFT framework which will be identi- 
fied soon, a = 1 — a denotes the "hole" density ma- 
trix. The involving matrix products are defined as usual; 
e.g., [C^tr] MV = E a =£,fl£ ! /Cii'(<vKv- Straight- 
forwardly, the current can be expressed in terms of the 
RSPDM as 



J(i) = eRe {Tt C^ ^(t) - C^ +) a(t) } 



(9) 



In arbitrary state basis, derivation is the same as above. 
The difference lies only at the expression of A^tJ , which 
in a non-eigenstate representation is given formally as 
Aa/j, = ^2i V C a iL V a v — ^ vv ^ m C atlu ,{^e m )D u , m D my a v . 
Here e m is the eigen-energy of eigenstate |to), and D is 
the transformation matrix from the non-eigestate rep- 
resentation to the eigenstate one. Obviously, with this 
identification, the resultant master equation and current 
formula are the same as Eqs. © and @, only replacing 
the matrices by C (± \ 

As an illustrative application of Eqs. (JSJ) and 10, we 
consider the simple non-interacting multi-level model 
studied in Ref.|8|. In the non-equilibrium stationary state, 
a(t — > oo) is diagonal in the eigenstate basis, thus [h, a] = 
0. As a consequence, the stationary state solution is de- 



termined by clfj (e^oy, 

ing to the well-known result [8|, c MM = [TL(e^)nL(e^j 
"R^Ffl^l/PM^) + T R (e^,)]. In particular, in the 
special case of equilibrium, tT MM reduces to the Fermi- 
Dirac function. Substituting er MM into Eq. ©, the well- 
known resonant tunnel current is obtained. 

(ii) Time Dependent System Hamiltonian: In this 
case, the RSPDM can be introduced in a similar 
manner. Consider, for example, Tr[a^j4[/ p(t)] — 

Y, a u>$ldt'C£ v S,t')°»'^',t) = [C<-VW. Here, 
t) = T!r{at[£(i,t')<v]/9(t)}, which can be solved 



- r {+) ( 



,), lead- 



via dfo-v'^t' \t) 



-i[h(t')a(t' ,t)]v>fj_, with the ini- 



tial condition ov^ (t, t) = Tr [a J, a v /p(t)]. Similarly, we 
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have Tr[Ai +) al p(t)] = £ Ql/ f* dt'C^, (t, t) = 

[C^a} u „. Here, a u> {t\t) = Tr{[a(*, i')<v]a+ p(t)}, sat- 
isfying an equation of the same form as a v < M (i' , i), but 
with initial condition d v i^(t,t) = 5 V *^ — <J u i^(t). As a 
result, in the time-dependent case, the resultant master 
equation and transport current can also be expressed as 
Eqs. © and J§J|, only keeping in mind that the matrices 
product needs not only the inner-state summation, but 
also the "inner-time" integration. 

Now we extend the above RSPDM formalism, i.e., 
Eqs.© and 0, to interacting systems. Within the 
TDDFT framework 01 > this can be straightforwardly 
done by replacing the single particle Hamiltonian by the 
Fock matrix 

hmn(t) = h° mn (t) + + <Tij(t)V mnij . (10) 

In first-principles calculation the state basis 
is usually chosen as the local atomic orbitals, 
{<fi m (r),m — 1,2, •••}. Here h°(t) is the non- 
interacting Hamiltonian which can be in general time- 
dependent; V mn ij is the two-electron Coulomb integral, 
V mnij = f dr J dr'r m (r)Mr)jdF]^(r')M^ and 
CW = ! dvcp* m {v)v^[n](v,t)4> n {v), with v™[n](r,t) 
the exchange-correlation potential, which is defined by 
the functional derivative of the the exchange-correlation 



functional A xc . In practice, especially in the time- 
dependent case, the unknown functional A xc can be 
approximated by the energy functional _E XC , obtained 
in the Kohn-Sham theory and further with the local 
density approximation (LDA). Notice that the density 
function n(r,t) appeared in the Fock operator is related 
to the RSPDM via n(r,t) = £ m „ <f> m (r)a mn (t)<f>l{r). 
Thus, Eqs. P- 1(10)1 constitute a closed form of TDDFT 
approach for the first-principles study of quantum trans- 
port, which is currently an intensive research subject 

To summarize, we have proposed a compact transport 
formalism from the perspective of quantum open sys- 
tems. The new formulation is constructed in terms of an 
improved reduced density matrix approach at the SCBA 
level, which is shown to be accurate enough in practice. 
Based on the established density matrix formalism, we 
also developed a new TDDFT scheme for first-principles 
study of transport through complex large-scale systems. 
Systematic applications and numerical implementations 
are in progress and will be published elsewhere. 
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